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We theoretically study the Shapiro steps in a hetero- Josephson junction made of a single-gap 
superconductor and a two-gap one. We find that an anomalous dc Josephson current is induced 
by tuning the frequency of an applied microwave to the Josephson-Leggett mode frequency, which 
creates an extra step structure in the I-V characteristics besides the conventional Shapiro steps. 
The step heights at the resonance voltages exhibit an alternate structure of a large and small value 
reflecting the gap symmetry of the two-gap superconductor. In the zbs-wave case in which the two 
gaps have opposite signs in the two-gap superconductor the steps with odd index are enhanced, 
whereas in the s-wave case the ones with even index have larger values. The existence of the 
fractional Shapiro steps is also predicted. 
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Since the discovery of iron-pnictide superconduc- 
tors^^^, the symmetry of the superconducting gaps has 
attracted a great interest. The spin- fluctuation mech- 
anism based on nesting between disconnected multiple 
Fermi surfaces predicts the ±s-wave symmetry It has 
been reported that the scanning tunneling microscope ex- 
hibits a characteristic field dependence expected from the 
dzs-wave symmetry^, and also the scanning SQUID de- 
tects half-fiuxon dynamics consistent with the ±5- wave 
gap^. However, the consensus on the pairing symmetry 
of the iron-based superconductors have not yet been at- 
tained^. Hence, a more definitive experimental probe is 
now in great demand. 

To clarify the pairing symmetry and associated new 
physics, various studies have been done for the Joseph- 
son effects in the iron-based superconductors both ex- 
perimentally^— and theoretically The present au- 
thors— have shown theoretically that a superconducting- 
insulator-superconducting (SIS) Josephson junction 
made of single-gap and multi-gap superconductors re- 
veals Josephson effects depending sensitively on the gap 
symmetry. In such a hetero Josephson junction the 
Cooper-pairs can transfer via multiple tunneling channels 
between the superconducting electrodes, because one of 
the electrodes has multiple conduction bands available 
for the tunneling. In the case of the ±5-wave symmetry 
the two tunneling channels, i.e., "0-junction" and "tt- 
junction" are formed. Therefore, one expects unusual 
Josephson effects different from the conventional one. 
Furthermore, a new phase oscillation mode called the 
"Josephson-Leggett (JL) mode" ^^'^^ exists in the low fre- 
quency region besides the Josephson plasma. This mode 
is expected to affect the AC Josephson effect. In this 
paper, we investigate effects of the JL mode on Shapiro 
steps in such an SIS hetero Josephson junction under 
microwave irradiation. We show that resonance between 



the JL mode and an applied microwave occurs in the 
phase running states when the microwave frequency is 
tuned to the frequency of the JL mode. The maximum 
DC Josephson currents induced at the resonant voltages 
exhibit remarkable behavior reflecting the gap symmetry. 
Furthermore, the existence of fractional Shapiro steps is 
predicted. 

Consider an SIS hetero Josephson junction shown 
schematically in Fig.[Tl We assume that the electrodes 
1 and 2 in this system are, respectively, single- and two- 
gap superconductors and the insulating barrier between 
the two electrodes has width d and dielectric constant e. 
To observe the Shapiro steps we irradiate a microwave 
generating the AC voltage l^wcosl^mw^ at the junction 
site as shown in Fig.[Tl 

Assuming that the phase difference is uniform along 
the in-plane direction (|| x?/-plane), one can derive the 
effective Lagrangian density in the present hetero Joseph- 
son junction system under no external magnetic fieldi^ 
as follows 

Sttw' ^ B-Kiii ^ e* 

i=l i=l 

+ ^cos^+^(^|i)2, (1) 

where 0^'^^ and if) are, respectively, the gauge- invariant 
phase difference and the relative phase of the supercon- 
ducting gaps having the phases Lp^^^ and (p^'^^ in the elec- 
trode 2 defined as 

with Lp^^^ being the phase of the electrode 1. The quan- 
tity A21 in 0^^^ is given as A21 — J^^f,^A^{z)dz, 
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FIG. 1: Schematic figure for a hetero Josephson junction 
made of single-gap and two-gap superconductors. A mi- 
crowave, whose amplitude and frequency are, respectively, 
Vmw and Qmw, is applied on the voltage-biased Josephson 
junction. 



A^{z) being the z component of the vector potential. 
The first and the second terms in Eq. ([T]) represent the en- 
ergy originating from the finite charge compressibility in- 
side the superconducting electrodes where the compress- 
ible area is confined around the interface whose thick- 
ness is characterized by the charge screening length ji' 
(for the electrode 1) or jii (for the electrode 2). We de- 
fine and qi as, respectively, = {fi/e*)dt(p^ + A\ and 
qi = {h/e*)dtip^^ + with e* = 2e and being the 
scalar potential in the electrode ^{= 1,2). The third 
term in Eq. ([T]) is the Josephson coupling energy in the 
two tunneling channels with critical current values ji and 
jf2. The forth term in Eq. ([T]) is the inter-band Josephson 
coupling energy, which originates from the Cooper-pair 
transfer between the two bands in the electrode 2. The 
coupling constant Jin can take both signs + and — . In the 
case of the ±5- wave (5- wave) we have Jin < (Jin > 0). 
Note that the inter-band coupling energy is minimum at 

= TT if Jin < 0. The fifth term in Eq. ([T]) is the energy of 
the electric field. The electric field inside the insulating 
layer is given as £^21 = —c~^dtA2i — d~^{A2 — A^). 

Let us study the dynamics of 9^'^\ From the effec- 
tive Lagrangian density ([1]) it follows the Josephson re- 
lation— 



where a' = j oti = ejii/d^ = a^^ -h cx^^ and 
A = 1 + a' + a. When the voltage appearing in the 
junction is given as (i^^fi — ^ + cos l^mw^, we obtain 



0(2) (t) = 



q(1) 

q(2) 
^0 



A/(i) 



Oil 
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(2a) 
(2b) 



us'mrtmwt, ^0 = e*Vo/h^ and u = 



where f{t) = Qot 

6^^^ = (tt) for Jin > 
(< 0). The quantity ip{t) describing the relative phase 



/HQrnw The phases ^q"^^ and ^q^^ are constants 
satisfying the relation, iIjq = Oq^-' 
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FIG. 2: DC components of /s(t) vs. Qo(= e^'Vo/h). The 
applied microwave has the frequency Qmw = cjjl (i.e., ^ = 1) 
and the amplitude u{= e*Vmw/^^^mw) = 1- The junction 
parameters are chosen as ^2/^1 = 0.93, a^lotx — 1.3, and 
ax = 0.1. (a) /3rit = 10"^ and (b) /3rit = 10"^ In (b), the 
steps labeled by even (odd) numbers are larger in the s-wave 
( ± s-wave) case. 



fluctuation, = ilj(^t) — ipQ. In this paper is assumed 
to be small. Since the total Josephson current in the 
present system is given as 



ls{t) = ^j,sm[0^'\t)], 



(3) 



2=1 



we have a simple expression, /g = jc('0o) sin[6>Q^^ + A/(t)], 

for tlj{t) = with jc('0o) = ji + j2cos?/^o- Then, us- 
ing the standard technique one finds that the Shapiro 
steps appear at the voltages which are the same as 
in the conventional Josephson junctions, i.e., Vq = 
— nA~^(fi/e*)l^niw, n being an integer. We note that 
the step height \jc{il^o)Jn{^u)\ depends on the relative 
phase, that is, the value of the critical current at each step 
is enhanced (suppressed) for the s-wave (zbs-wave) case. 
The critical current depending on the gap-symmetry also 
leads to the anomaly of the Reidel peak, which can be 
used for the identification of the ±5- wave gap as proposed 
by Inotani and Ohashi^^. 

Next suppose that ?/^(t) is temporally fluctuated. To 
describe the time evolution of we utilize the equation of 
motion obtained from Eq. ([1]). From the Euler-Lagrange 
equations with respect to (p^ and one can derive the 
equation as — , 
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i) 
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-cji sm 



^(1) 



;sin(9(^) - sgn(Jin)2a;f^sinV^, (4) 



where Ui = y^47re*dji/eh, cjin = Y^47re*(i| Jin|/efi, and 
^ = (<^i — 0L2)/{ai + a2). We note that Eq. (j4]) describes 
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FIG. 3: DC components of /s(t) vs. Qo(= e*Vo/h). The fre- 
quency and the amphtude of an apphed microwave, respec- 
tively, Qmw = ^cjjL {q / 1) and u{= e*Vmw/^^mw) = 1- The 
detuning P^it = 10~^. {a) q = 2 (b) ^ = 3. 



the effect that the charge imbalance between the two 
bands in the electrode 2 (the left-hand-side of the equa- 
tion) is induced by the Josephson current between the 
two electrodes (the right-hand-side). The charge imbal- 
ance in this system is characterized by the parameter ^. 
Substituting Eqs. (j2aj) and (j2b]) into Eq. (j4]), we find the 
linearized equation for small , which is similar to the 
one for a forced harmonic oscillator, 



5tV = -['^j'l + Gjim + Fcit) + Flit), 



(5) 



where 



Fc{t) -^^^-^^nl^^usinnrnwt, 



{cvj+ojl)sm[ei'^ +Af{t)], 



with u = [uJi— sgn(Jin)cj2]/(^i + ^2) Gj{t) = 

Yl'i=i{(^i^V^) cos[6>^'^ + A/(t)]. The external force terms 
Fc and Fj comes from, respectively, the charge imbalance 
and the Cooper pair tunnehng. The angular frequency 
cjjL is given as cjjl = + Q^2'^in, which is equal to 
the mass of the JL mode—. We assume |Jin| ^ ji, j2 
throughout this paper. In this case, noting uf^ 
one can safely neglect Gj in Eq. (|5]). Then, under this 
assumption we have a solution of Eq. (|5|) as follows^^, 

^{t) = Asmnrn^t^"^Bks'm{nkt^e^^^), (6) 
where the coefficients A and Bj^ are given as 



Bk = -1^ 



2 

jl+j2 ^fe(A^) 
2 Jin 



with Clk = Allo+^^mw and Jk{^u) being the Bessel func- 
tion of the kih. order. It should be noted that 7/^ has poles 
at 1]^^ = cjjL and Vt\ = cJjl, that is, i/j is resonantly en- 
hanced at these frequencies. These resonant conditions 
can be satisfied by tuning the frequency of the applied 
microwave. In a region far from the resonance we have 
ip 0. In the region we have conventional Shapiro steps 
as discussed above. It should be also noted that Eq. (j5j) 
has a characteristic oscillating solution corresponding to 
the JL mode, which is dropped in Eq. (|6|) because the res- 
onance effects are predominant. An enhancement of the 
DC tunneling current by this collective mode is pointed 
out without microwave irradiation as discussed in Ref. 

Let us now study the resonance effects due to the 
microwave irradiation on the I-V characteristics in the 
present system. It is convenient to use the dimensionless 
parameters, p and defined as 



qoJjL- 



The frequency is then expressed as = {p-\- k)qujji,. 
Close to the resonance we write the coefficients A and 
5fe as A = -^[(ai + a2)/2]{u/f3rit) and Bk = + 
j2)/2Jin][J/c(Ai^)//3rit], where P^it is regarded as the de- 
tuning parameter. Let us first examine the case of g' = 1 
and k = k± = —p ±1. In this case the resonance takes 
place at flmw = ^jl and ftk = i^jL, and thus Eq. ([6]) is 
approximated as ~ Asinc^jL^ + B^^ sm{ujji,t + ^q^^) + 

Bk_ sin(— cjjL^ + ^0^^)- Then, substitution of Eqs. ([2a|) 
and (|2bp into Eq. ([3|) and use of the Fourier-Bessel expan- 
sion for Eq. ([3|) leads to a condition that a DC compo- 
nent emerges in Is{t), i.e., AQq + (n + — £_)cjjl = 0, 
where n,i± G Z. Hence, the voltages at which the 
Shapiro steps appear can be determined by the relation, 

p{q = 1) = -n-^+ +^_(g Z). 

Let us present the calculated results at g = 1. We 
choose the parameter values as u = 1, j2/ji = 0.93, 
OL2loL\ = 1.3, and a — 0.1. Figure [2] shows the Shapiro 
steps in the cases of (a) /3rit = 10~^ and (b) /3rit = 10~^. 
As seen in this figure, the voltage dependence of the 
step height drastically changes as one approaches the 
resonance condition, i.e., from a monotonic decrease to 
oscillating one. Furthermore, one notices in Fig.[2fb) 
that in the oscillating region a larger DC current ap- 
pears at the voltages corresponding to even (odd) p in 
the case of the s-wave (±5- wave), that is, the Shapiro 
steps generated by the resonance with the JL mode 
clearly depends on the pairing symmetry. To under- 
stand the origin of this remarkable feature in the Shapiro 
steps we investigate the induced DC current for large 
p. It is easy to show that ip ~ Asinc^jL^ for |p| ^ 1 
— . Then, the step height in this limit is obtained as 
^s,DC = |ji-^p(^i) + j2CosV^o-^p(2:2)|, where xx = K.u ^ 
{pL\A) I {oL\ + 0L2) and X2 = Aix — (a2 A) / [ol\ + 0^2). Since 
A is large near the resonance frequency, one may have an 
approximate relation, x\ o:^ —X2- In this case, employing 
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the relation J^(— x) = (— l)^Jn(x) (n G Z and x G M), 
we find that 

^s,DC - b'l + (-lFj2 COS V^ol • 

The coefficient \ji + (— l)^j2 cos?/^o| takes alternating val- 
ues with respect to p, that is, ji + (— l)^j2 in the 5- wave 
case {ipo = 0), whereas ji + (— l)^^"^j2 in the ±5- wave 
case {i/jQ = TT ), which leads to the fact that a larger DC 
current appears at even (odd) p in the 5-wave (±5-wave) 
case. This remarkable feature of the Shapiro steps caused 
by the resonance with the JL mode is quite clear-cut for 
the determination of the pairing symmetry in multi-gap 
superconductors. 

Finally, we examine the case of q > 1 {q G Z). The 
resonance in this case arises from the poles in the second 
term on the right hand side of Eq. (|6]), i.e., flk = ^jl, 
which brings about the correction to 6^'^\t) as '0 ^ 
sin(u;jLt + ^q"*^^). Then, after similar mathematical 
manipulation to that in the q = 1 case one finds the re- 
lation. Alio + ^^mw + ^/c^JL = with n, ^k ^ Z, for the 
DC current to be induced. From this result it follows, 

p{Vtk ^JL, g > 1) = -n - — = -/c + -. 

Q q 

We note that p is generally not an integer, which indicates 
that the Shapiro steps appear at fractional voltages, that 
is, the resonance with the JL mode induces the fractional 
Shapiro steps, as well. Figure [3] shows the DC current. 



i.e., the step heights at the resonance voltages. The mi- 
crowave amplitude and the junction parameters are the 
same as in Fig.O Figures [3l^a) and (b) show the cases 
of q = 2 (i.e., p = m ^ ^ (m G Z)) and q = 3 (i.e., 
p = m + I (m G Z) ), respectively. The numerical results 
given in these figures proves the existence of the fractional 
Shapiro steps. Although fractional order steps have been 
studied in Josephson arrays and highly transmissive 
junctions with conventional one-gap superconductors, 
the origin in the present system is an intrinsic feature of 
multi-gap superconductors, the JL modes. 

In conclusion, we revealed that an external microwave 
applied to the hetero Josephson junction resonantly ex- 
cites the JL mode and induces DC currents. The volt- 
age dependence of the Shapiro steps is sensitive to the 
inter-band sign change, i.e., the pairing symmetry of the 
multi-based superconductor. It is also predicted that the 
fractional Shapiro steps appear. We suggest that the tun- 
ing of the microwave frequency to the JL mode one bring 
about not only unconventional types of Shapiro steps but 
also a clear probe to the gap symmetry of a multi-gap su- 
perconductor. 
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